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Abstract 

The gluon spin coupling to a Gaussian correlated background gauge 
field induces an effective tachyonic gluon mass. It is momentum de- 
pendent and vanishes in the UV only like 1/p 2 . In the IR, we obtain 
stabilization through a positive mconfiP 2 ) related to confinement. Re- 
cently a purely phenomenological tachyonic gluon mass was used to 
explain the linear rise in the qq static potential at small distances and 
also some long standing discrepancies found in QCD sum rules. We 
show that the stochastic vacuum model of QCD predicts a gluon mass 
with the desired properties. 



Recently it was argued |1|, that an effective negative gluon (mass) 2 
produces the linear rise of the static quark-antiquark potential observed in 
lattice calculations |§ at rather short distances below the nonperturbative 
IR scale of QCD. This can also be related to an unconventional Q~ 2 piece 
in the running QCD coupling a s {Q 2 ) j|, |l| indicated again by lattice data 
on the gluon condensate and the 3-gluon vertex ||. Later in ref. 0, it was 
shown that a linear term in the potential for small distances appears in the 
evaluation of the Wilson loop integral in a Gaussian correlated background 
gauge field if the paramagnetic (spin related) gluon coupling is taken into 
account. A gluon propagator with a nonperturbative tachyonic mass in the 
evaluation of the operator product expansion for QCD sum rules also solves 
old discrepancies in the pion and scalar gluonium channels 0. 

In ref. 0, the view was taken that a gluon mass is not just a convenient 
way to parametrize nonperturbative infrared (IR) effects, but that a tachy- 
onic gluon mass also enters into the basic short distance operator product 
expansions. In this letter, we propose less radically that a tachyonic gluon 
mass does appear at momenta above the QCD scale, although at very high 
momenta/small distances it goes to zero. Thus, in the real ultraviolet regime 
(UV), the basic perturbative operator product expansion should be valid|] 
Of course, if one discusses an operator product expansion in an effective the- 
ory including a tachyonic gluon mass in which very high momenta are already 
integrated out, this mass shows up in the expansion. However, at very small 
distances, the expansion loses its basis. 

In the IR, a tachyonic gluon mass would cause divergences. Thus there 
should be an overcompensating positive effective gluon (mass) 2 . In the case 
of the hot electroweak theory close to the symmetric phase, we have argued 
earlier || for the following picture: a tachyonic gauge boson mass related to 
spin-spin forces is effective at intermediate momenta and a positive (mass) 2 
related to the area law of confinement is dominant in the gauge boson propa- 
gator in the IR. Our equations in ref. |5| were written for a general dimension 
d and gauge group SU(N), but we had restricted our discussion to the case 
d — 3 and N = 2, which is the relevant one for the discussion of the elec- 
troweak phase transition. In this letter, we analyze the QCD case d = 4 and 
N = 3. 



■"Concerning a l/Q 2 behavior consider the following simple function: (aQ 4 + (3Q 2 ) 
If 3> a is large, it behaves like -gkj for Q 2 < (3 /a, i.e. for rather large Q 2 \ 



Our actual calculations use the "stochastic vacuum model" JTIJ for the 
QCD vacuum. In brief, this model provides the area law of confinement 
in a very simple and intuitive way. It describes the QCD vacuum by a 
stochastic process and assumes that only Gaussian correlations are present 
in the cumulant expansion. The nonperturbative part of the gauge invariant 
field strength correlator is given by (after contracting the Lorentz indices): 



((9 2 FZ„(x',x )F» u (x,x ))) = 6 ab (g 2 F 2 )D (^-^) • (1) 



Here (g 2 F 2 ) = (g 2 F^ l/ F^ u ) is the usual "local" gluon condensate, and a 
denotes the correlation length. The tensor F iiv {x 1 xq) is the ordinary field 
strength tensor F^x) parallel-transported to a fixed common reference point 
Xq pi and a Euclidean metric is used. A typical form of D used in lattice 
evaluations [111] is exp(— \x' — x\/d). 



The inverse gluon propagator (in the Feynman background gauge) in a 
background gauge field reads 

K% = [-D 2 <V + m 2 l c <V + 2igF liV ] *" . (2) 

Here we have allowed for the presence of a gluon mass for future convenience. 
Linearized gauge field excitations a^(x) interact with the background field via 
the Lagrangian ^a^K^al. We have to distinguish the "diamagnetic" cou- 
pling via the ordinary minimal-substitution term D 2 5^ u and a non-minimal 
"paramagnetic" coupling via the 2igF fll/ term. In || we used this interaction 
Lagrangian in order to compute the "mass operator" E of the gluon (vacuum 
polarization) to lowest nontrivial order in (g 2 F 2 ). It is defined such that 

G= [ G (o)-i_ s] -i ( 3) 

where G^ and G are the free and the dressed propagator of the gluon, 
respectively. The latter describes the propagation of a perturbative gluon in 
the background of a Gaussian-correlated random Yang-Mills field. Taking 
xq = (x + x')/2 as the reference point, the result for E is diagonal in the 
Lorentz and the color indices: 

E(x, x')% = S{z 2 = {x- x') 2 )5, u 5 ab . (4) 

The function S(z 2 ) = Sp(z 2 ) +S^(z 2 ) receives contributions from two rather 
different physical mechanisms: Sp results from the spin-related paramagnetic 



interactions via the F^ term alone, while Sa stems from diagrams involv- 
ing the diamagnetic coupling and coincides with the result that we would 
obtain for spin-0 bosons. Fourier-transforming with respect to z, we obtain 
S(p 2 , m 2 ) = S F + S A withfl 
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(5) 



and similarly for Sa{p 2 , m 2 ). 

In view of eq. (0) it is natural to interpret fn| ac h(p 2 ) — —S(p 2 ,m 2 = 
0) as the effective, momentum-dependent mass of a gluon propagating in 
a stochastic gauge field background. One finds that S{p 2 , 0) > so that 



m 



tach 



< 0, i. e. the effective mass is tachyonic. Its value at p 



m tach(0)> can be written down in closed formP]: 

3 



0) m tach 



- mU = S(p 2 = Q,m 2 = Q)= a 2 {g 2 F 2 ){l + 5} dw D(w 



(6) 
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Here the "1" inside the curly brackets on the RHS of eq. (f|) comes from the 
paramagnetic Sp while 
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-0.492 



5 = 1Z 

12 3 

is the contribution of the diamagnetic piece Sa- We observe that the tachy- 
onic nature of the induced mass is a consequence of the paramagnetic spin 
interaction which dominates over the diamagnetic one. The latter generates 
a (mass) 2 which is numerically smaller and of opposite sign. 

It is quite remarkable that, up to a constant factor, the RHS of eq. (M) 



equals precisely the well-known |10[ result of the stochastic vacuum model for 
the string tension (with the color sources in the fundamental representation 
of 577(3)): 

(7) 



7T f°° 

Cfund = -—a 2 {g 2 F 2 ) / dwD{w). 
144 Jo 



tSee eq. (A.34) of 
*Cf. eq. (A.38) of [ 



where w = z 2 /a 2 
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Figure 1: SF{m 2 or£ = 0) (solid line), SF{rn 2 on{ 7^ 0) (dashed-dotted line) and 
m conf (dashed line) as function of p 2 . 



Thus the gluon mass can be expressed directly in terms of the string tension: 



2 18 „ 5 . 

™tach = — (hi2-— a fU nd 
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2.18(Tf un d. 



(8) 



It is a universal prediction of the stochastic vacuum model in the sense that 
no separate knowledge of the gluon condensate, the correlation length and 
the shape of the structure function D is needed. In fact, in eq. (|J) we 
can directly use the experimental value for the string tension. From the 
charmonium spectrum one obtains |D| crf un d ~ (4:30MeV) 2 which leads to 
m 2 ach « -0A0(GeV) 2 or |m tach | « 635MeV. 

This is indeed a tachyonic mass of the order of magnitude postulated in 
ref. || on phenomenological grounds. Since to lowest order 
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with an l/p A correction to the free propagator, we see that in the static 
qq potential V(r) the tachyonic gluon mass gives rise to a term that grows 
linearly with r. However, in contrast to the asymptotic (r — > 00) regime, in 
which V(r) ~ <J{ un dr dominates the Coulomb term, we are at much shorter 
distances here where the linear term is only a small correction to the 1/r law. 
In fig. 1, we present SF(p 2 ,m 2 = 0) for D(w) = exp(— \z\/a) = exp(— \w\i) 



and a = 0.22fm [O] as a function of p 2 . It vanishes for large p 2 rather slowly 



as 

s F ^ -, §££>. (10) 

Eqs. (f|), (|]) and (|10|) are our main results. The "paramagnetic domi- 
nance" and the resulting tachyonic sign of the gluon (mass) 2 are related to 
the well known instability of the Savvidy vacuum ||14| . The latter models 
the true QCD vacuum in terms of an external gauge field with a constant 
field strength. While already lower in energy than the perturbative vacuum, 
it is unstable and decays towards the true ground state precisely because 
the destabilizing paramagnetic interaction of the vacuum fluctuation over- 
rides the diamagnetic one. The paramagnetic dominance encountered above 
(and discussed in detail in ref. 0) is a remnant of this phenomenon for 
non-constant backgrounds. 

The self-energy £ can also be used in order to understand why, first of 
all, the perturbative vacuum is unstable. Inserting E, which is proportional 
to (g 2 F 2 ), into a perturbative gluon loop, we obtain the (F 2 ) term of the 
1-loop effective potential for the condensate, V e s({F 2 )). This term turns out 
to be negative, i. e. (g 2 F 2 ) = is not the minimum of V e s, and a non-zero 
condensate tends to form. But, in contrast to the d = 3 case, for d = 4, the 
term that is linear in (F 2 ) is logarithmically divergent in the UV. 

In ref. @ , we also calculated the standard perturbative 1-loop gluon self- 
energy II and we introduced the F^ v correlator (|l|) only in a second step (fig. 
2). This should be an equivalent procedure and indeed, in three dimensions, 
this can be easily checked. For d = 4, the perturbative vacuum polarization 
II(g 2 ) (fig. 2) is divergent. It has to be renormalized at some Euclidean 
q 2 = fi 2 in the perturbative range /i 2 > l/a 2 .| Taking only the contribution 
from the non-minimal gauge boson coupling, one obtains (see eq. (2.20) of 
ref. i) 

n ^ 2 '^ ) = (4^ l0 V (U) 

which is to be inserted into 

V« = \(9 2 F 2 ) J ^Wr( P 2 ,^)D cS (p 2 ). (12) 

§ According to ref. p2| the stochastic vacuum model describes the QCD vacuum at a 
renormalization scale fi not much above 1/ a. 



Figure 2: First order contribution of the correlator (double wiggled) to the 
gluon loop 

This is indeed negative if /x 2 > l/a 2 , the typical scale set by D. 
The 1-loop expression resulting from the £ approach, 

Va ex J ^i(log(p 2 - MP 2 )) - logp 2 }, (13) 

with an appropriate constant already subtracted, is still logarithmically di- 
vergent. Obviously the 1/p 2 decay of S(p 2 ) is too slow to make £ completely 
irrelevant in the UV. At first sight, this seems to be an UV-effect related to 
the gluon mass similar to the one discussed in ref. @. However, our divergent 
term is proportional to the stochastic average of the classical action \F^ u F^ y 
and can be removed by the same renormalization as in perturbation theory. 
This is an important difference as compared to ref. [J2J. Note also that in 
contrast to the case with a constant background field, we are not plagued by 
IR singularities any longer. 

To really obtain an instability at (F 2 ) = 0, the positive tree term has 
to be dominated by the quantum corrections. The 1-loop potential could be 
further improved by including Wilson type renormalization effects [|HJ: the 
effective g 2 (k 2 ) is increasing in the infrared direction k 2 — > 0. 

Finally, let us discuss the gluon propagator and the stochastically aver- 
aged 1-loop effective potential 

V eS ((F 2 )) = i((Tr log(J0» (14) 



beyond the linear approximation in (F 2 ). The effective potential can be 
represented in terms of a world-line path integral |J over closed paths y(r). 
Its integrand contains the pivotal factor 

{(tT^Pexpig J drfaAM + 2F(y)))). 

Apart from the paramagnetic F-term discussed above, this is precisely the 
stochastic average of the Wilson loop operator which signals confinement. 
For large loops, it behaves as exp(— aA[y\T) where a = o"f un d oc (g 2 F 2 ) is the 
string tension ([?D and A[y] denotes the area of the minimal surface bordered 
by y(r). This is one of the most important results of the stochastic vacuum 
model. Thus the contribution to V e g which originates from the "diamagnetic" 
interaction of very large loops is proportional to 

1 r°° rIT r r 1 

V ^ = -2J -Y T 'j DyeM ~L dT y 2 /^){^v{-<yA[y]Tf{A[y}T))-l}. 

(15) 
Here we use appropriately rescaled dimensionless variables. The path integra- 
tion is subject to the conditions y(0) = y(l) and J^ y(r)dT = 0. In principle, 
the cut-off function / with /(oo) = 1 and /(0) = could be calculated from 
the model; it serves the purpose of cutting out the contribution of small loops 
so that the perturbative expression at small distances is not changed. Using 
the identity (4irT)~ 2 = J -S4% exp(— Tp 2 ), the confinement-related part of 

Kfr assumes the form[] 

Kr«^ 2 »= I-^tAp 2 ) ( i6 ) 



and 



(27T) 

I(p 2 ) = -4vr £J du>F{u) f ^ eM~P 2 T)G (£) (17) 

F(u) = J Dy exp \iA[y]u - £ dry 2 /^ (18) 

/■oo _ _ r _ _ -, 

G(x)=2 dAcos{Ax) exp{-aAf(A)} - 1 . (19) 

Jo L J 



where 

'"+ 00 , „, s f°° dT , ~ s „ (uj 



with 



"One of the authors (M. G. Schmidt) thanks M. Laine for a helpful discussion on that 
point. 
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Expression (fHJ) can be written in the form {ln(p 2 + m 2 onf (p 2 )) — logp 2 }. 
Eq. (|16|) suggests that this defines an effective IR mass fnl on {(p 2 ) due to 
confinement ||. Note that the counterterm cancels in this relation. 

The evaluation of (|TTD, (|18|) and (p~9f) requires numerical studies, in par- 



ticular of the most interesting function F(u) [ 16| . The main scale in this 
problem is the string tension a; thus we expect m 2 orii (j) 2 = 0) ~ a. The 
shape of the cut-off function f(A) determines the profile in p 2 . At p 2 = 0, 
m 2 onf should dominate the negative —Sf{j> 2 = 0) in order to obtain a reason- 
able IR behavior (though interferences between confinement and spin effects 
are not treated carefully in the approach presented above). 

The confinement mass m 2 oni {p 2 ) should be also included in our calculation 
of Sp- In a first approximation, one could substitute it as an m? in formula 
(ID (of course, it is again not an ultraviolet mass requiring regularization). 
This lowers Sp and the result depends a lot on the balance between m 2 onf and 
—Sp at small p 2 . In fig. 1, we have included a plausible function m 2 onf (p 2 ) to 
demonstrate the effect. We expect m 2 on {{p 2 ) to dominant in the IR, i. e. for 
p 2 < 1/a 2 . The range of Sf(p 2 ) extends to higher p 2 . It vanishes only like 

1/p 2 . Thus —Sf{p 2 ) mimics a tachyonic UV mass. 
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